Diffusion of Neutrinos in Proto-Neutron Star Matter with Quarks 
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Abstract 

Neutrino opacities important in the evolution of a proto-neutron star containing quark 
matter are studied. The results for pure quark matter are compared with limiting 
expressions previously derived, and are generalized to the temperatures, neutrino de- 
generacies and lepton contents encountered in a proto-neutron star's evolution. We 
find that the appearance of quarks in baryonic matter drastically reduces the neu- 
trino opacity for a given entropy, the reduction being sensitive to the thermodynamic 
conditions in the mixed quark-hadron phase. 

PACS: 97.60.Jd, 21.65.+f, 13.15.+g, 26.60.+C 



A general picture of the early evolution of a proto-neutron star (PNS) is becoming well established 
E H . Neutrinos are produced iu .ar g e q ua„ titl e S by eiectrou capture a 8 the proper star eo»a p8 e 8 , 

a ■ 

but most are temporarily prevented from escaping because their mean free paths are considerably 
smaller than the radius of the star. During this trapped-neutrino era, the entropy per baryon s is 
about 1 through most of the star and the total number of leptons per baryon Yl = Y e + Y Ue ~ 0.4. The 
neutrinos trapped in the core strongly inhibit the appearance of exotic matter, whether in the form of 
hyperons, a Bose (pion or kaon) condensate or quarks, due to the large values of the electron chemical 
potential. As the star cools, the neutrino mean free path increases, and the neutrinos eventually leak 
out of the star, on a timescale of 20-60 s. During deleptonization, neutrino diffusion heats the matter to 
an approximately uniform entropy per baryon of 2. If the strongly interacting components consist only 
of nucleons, the maximum supportable mass increases. In the case that hyperons, a Bose condensate 



(pion, kaon) or quarks appear in the core as the neutrinos leave, the maximum mass decreases with 
decreasing leptonic content. Neutron stars which have masses above the maximum mass for completely 
deleptonized matter are thus metastable, and will collapse into a black hole during deleptonization. 
Alternatively, if the mass of the neutron star is sufficiently small, the star remains stable and cools 
within a minute or so to temperatures below 1 MeV as the neutrinos continue to carry energy away 
from the star. 

The way in which this picture is modified when the core of a PNS contains deconfined quark 
matter is only beginning to be investigated 0j7|-||. I n n ^ s seminal paper, Iwamoto [10] noted that 



the non-degenerate v mean free path in cold quark matter is about ten times larger than in nucleonic 
matter. We find that in PNS matter, in which quarks appear towards the end of deleptonization, 
similarly large enhancements persist even up to the largest relevant temperatures (~ 30 — 40 MeV ||), 
inasmuch as quarks remain largely degenerate. On this basis, it can be anticipated that the presence of 
quark matter increases the neutrino fluxes while simultaneously decreasing the deleptonization time, 
relative to matter without quarks. In work to be reported elsewhere, we explore the possibility that 
such a change might be detected from a Galactic supernova in current and planned neutrino detectors. 
This would have direct implications for the theoretical understanding of the high-density regime of 
QCD which is inaccessible to high energy Relativistic Heavy-Ion Collider experiments, and, currently, 
to lattice QCD calculations at finite baryon density. 

To perform detailed simulations of the neutrino signal from a PNS containing quark matter, as has 
already been done for matter containing nucleons, hyperons and/or a kaon condensate consistent 
calculations of neutino interactions in hot lepton-rich matter containing quarks are required. It is most 
likely that quarks exist in a mixed phase with hadrons [fF||i||TT| . Steiner, et. al. |§, recently showed 
that the temperature of an adiabat decreases as a function of density in a mixed phase of quarks and 
nucleons. Because v— cross sections usually scale with T 2 , this suggests that the presence of quark 
matter might influence the neutrino signal of a PNS with quarks. 

In this work, we calculate the diffusion coefficients of neutrinos in a mixed phase of hadrons and 
quarks for the temperatures, neutrino degeneracies, and lepton contents likely to be encountered in the 
evolution of a PNS with quarks. We demonstrate that the cross sections for scattering and absorption 



of neutrinos by nucleons, leptons, and quarks are reduced to two integrals, whose integrands are 
products of simple polynomials and thermal distribution functions. The limiting behaviors of the 
cross sections, for non-degenerate and degenerate neutrinos, respectively, are compared with previous 



calculations [T(J in the case of pure quark matter. For simulations of PNS evolution using the diffusion 
approximation, diffusion coefficients, which are energy weighted averages of neutrino cross sections, 
are required for matter in which quarks exist in a mixed phase with hadrons. We examine the relevant 
diffusion coefficients for two thermodynamic conditions especially germane to PNS evolution. The first 
situation is when neutrinos are trapped, s«l, and the total lepton content of the matter Yl = Y e + Y Ve 
(which measures the concentrations of the leptons per baryon) is approximately 0.4. We also consider 
the situation when neutrinos have mostly left the star and the matter has been diffusively heated 
(Y„ ~ 0,s ~ 2). We discuss the impact these results might have upon the evolution of a PNS which 
contains quarks in a mixed phase. 

For the v— energies of interest, the neutral and charged current v— interactions with the matter in 
a PNS are well described by a current- current Lagrangian jnj 

C = S (^(l)(l _ 75 )^ 3 (3)) (> P2 (2)(V - ^7 5 )VV(4)) + H.C. , (1) 

where V and A are the vector and axial-vector coupling constants (See Table 1) and Gf — 1.17 GeV~ 2 
is the Fermi weak coupling constant. The subscripts i — 1, 2, 3, and 4 on the four-momenta pi denote the 
incoming neutrino, the incoming lepton, baryon or quark, the outgoing neutrino (or electron), and the 
outgoing lepton, baryon or quark, respectively. The charged current reactions contribute to absorption 
of neutrinos by baryons or quarks and scattering of neutrinos with leptons of the same generation. 
The neutral current interactions contribute to scattering of neutrinos with leptons and baryons or 
quarks. The charged current contribution to neutrino-lepton scattering in the same generation can be 
transformed into the neutral current form, which modifies the constants V and A for that case. For 
completeness, the values of V and A for electron neutrinos are given in Table 1. The corresponding 
values for reactions with electron anti-neutrinos are obtained by the replacement A — > —A. 
From Fermi's golden rule, the cross section per unit volume (or inverse mean free path) is 

^(1-A)(1-/<)(2,)V< PI+P2 -P3- P1 ), (2) 



where the degeneracy factor g is 6 (3 colors x 2 spins) for reactions involving quarks while it is 
2 (2 spins) for baryons of a single species. The Fermi-Dirac distribution functions are denoted by 
fi — 1 + exp ( El ~^ ) , where E^ and /ij are the energy and chemical potential of particle i. The 
transition probability Wk, summed over the inital states and averaged over the final states, is 



W 



fi 



Or p 



(V + Af ^ ■ p 2 ) (p 3 ■ p 4 ) + (V - -4) 2 (pi ■ p 4 ) (P3 • P2) - (V 2 - A 2 ) (pi • pa) (P4 " P2)] • (3) 
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Utilizing <i 3 pj = p? dpidQi = piEi dEidVti and integrating over P 4 , Eq. fl2|) may be cast in the form 
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- = ^r T dE 2 T~ dE 3 S ^Nlpll (V + ^l) 2 4 + (V - A) 2 I b + (V 2 - A 2 ) I c 
V S2ir b Jm 2 Jo Ei L v 7 

where S — f 2 (1 — ^3) (1 — / 4 ), Mj is the mass of particle i, and 

Ia = J rf0 2 rf0 3 cZ04 5 3 (pi + p 2 - P3 - P4) (Pi ■ P2) (P3 • P 4 ) • 



(4) 



(5) 



The integrals I& and J c are defined similarly to I a , and all can be performed analytically. Explicitly, 
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(6) 



where 



2A = 2E 1 E 2 +p 2 1 +pl 2B = 2E 3 E 4 + pl+pl, 
P min = max(|pi - p 2 |, |p 3 -p 4 |) , Pmax = min(pi +p 2 ,p 3 + P 4 ) • 



(7) 



In the above expression, pi = p%. h and I c are defined to be the same as I a , but with appropriate 
replacements: 



h = Ia(P2^P4,E 2 ^-E 4 ), I c = I a (p 2 ^p 3 ,E 2 ^E 3 ). (8) 

In Eqs. © and (§), P 4 = Pi + E 2 - E 3 and |p 4 | = yJ(E 1 + E 2 - E 3 ) 2 - Mf . Eqs. (§ through © 
are the principal results of this work and enable us to compute, for arbitrary conditions of neutrino 
degeneracy and matter's temperature, the neutrino diffusion coefficients required in simulations of PNSs 
containing quark matter. We note that similar techniques were employed to calculate v— emissivities 
in cold catalyzed neutron stars in Refs. [|13"1 , |14| . 



In limiting cases when the neutrinos are either degenerate or non-degenerate, and the quarks, which 
are always degenerate in PNSs, are massless, simple analytical expressions for the cross section may be 
obtained by replacing momenta by Fermi momenta and energies by chemical potentials in the integrals 
I a , h, and I c . For the sake of comparing such limits with the general results obtained from Eq. (|]), 
we record various limiting forms obtained earlier in Refs. [|TU|,|i7) . 



(1) Scattering of degenerate neutrinos: The result is the same as for neutrino-electron scattering: 
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V 2 + A 2 ) (10 + x 2 ) +5(2VA)x 



(9) 



where x = min(Ei, /X2)/niax(Ei, /i2). Here, and in the following, we have removed the factor 1 — /1 
(l + e~( El ~^ 1 ^ T ^J from ||l0| , |i"5f to obtain transport mean free paths. 
(2) Scattering of non-degenerate neutrinos: The inverse scattering mean free path is 

a s G F E\& 



V 5tt 3 ' 

when it is additionally assumed that H2 is large compared to E\. 
(3) Absorption of degenerate neutrinos: 

2G 2 F fil 



(10) 
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5tt 3 [i 2 



(lO/i 2 + 5/i 4/ u 3 + jj 2 ^ {Ei - /ii) 2 + 7r 2 T 



(11) 



(4) Absorption of non- degenerate neutrinos: The general result is greatly simplified by additionally 
assuming that the quark chemical potentials are modified by perturbative gluon exchange: 

o A 16 



'-a c G' F pF 2 p F3 p F4 E{ + 7r 2 T 2 



(12) 



V 7T 

Representative cross sections from Eq. (|j) are compared with the limiting forms in Eqs. fl9|)-(|T2"l) 
in Fig. p]. Degenerate neutrinos are assumed in this figure to have fi u » T and non-degenerate 
neutrinos are assumed to have [i u « 0. In the regions where they were expected to be valid, namely 
E v jT » 1 for degenerate absorption and scattering, and also non-degenerate scattering, and E v w T 
for non-degenerate absorption, the limiting forms give adequate representations of the general results. 
However, significant deviations occur in the cases of non-degenerate absorption when E v 7^ T, and 
for degenerate absorption and non-degenerate scattering when E v << T. The deviation for non- 
degenerate absorption is due to the neglect of p v in the momentum conservation condition in Eq. (pj) in 



Ref. [|10j (which is appropriate for cold catalyzed stars, but not for hot matter in PNSs), which limits 
its applicability to the region E v « T. The other two deviations are simply due to the assumption in 
Ref. that E V »T. 

The weak interaction timscales of neutrinos are much smaller than the dynamical timescale of PNS 
evolution, which is on the order of seconds. Thus, until neutrinos enter the semi-transparent region, 
they remain close to thermal equilibrium in matter. Hence, neutrino propagation may be treated in the 
diffusion approximation with the differential equations for the flux of energy (F u ) and lepton number 
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and 



dr ^ ) A dr 

where A and <fi are general relativistic metric functions, rj = fi u /T and D2, -D3, and D4 are diffusion 
coefficients decomposed as 



D 2 = D^ + Dl% 



D 3 = D? 



-^3 1 



(13) 



The transport of \i and r neutrinos and anti-neutrinos are well approximated [|I|] by assuming that 
they contribute equally to and are represented by D v ». These coefficients are defined in terms of 
the energy dependent diffusion coefficient D P (E U ) by 



dx x n D p (E 1 )f(E 1 ) [1 - fiEj] 



(14) 



where x = Ei/T, and the superscipt p denotes either the electron neutrino, the anti-electron neutrino, 
or the /j, and r neutrinos and thier antineutrinos. In turn, the energy dependent diffusion coefficient is 
obtained directly from the cross sections per unit volume through 

1 
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E y + x E £ + (i-x) E % 

r=(p,L) V r=(p,H) r=(p,Q) V 



(15) 



where the (p, L), (p, H), and (p, Q) represent the sum over all the reactions of particle p with leptons, 
hadrons, or quarks, respectively. The factor (1 — /1) 1 ensures detailed balance, and x is the volume 
fraction of matter in the hadronic phase. 



Note that in the case of scattering, the Pauli blocking factor corresponding to the outgoing neutrino 
is omitted, since the neutrino distribution function is not always known a priori unless a full transport 
scheme is employed. It is possible, however, to devise a simplified scheme [[J in which the dependence 
on the neutrino distribution function is minimized. Such a scheme is valid only when scattering from 
light particles (either electrons or quarks) does not dominate the opacity. Our results below show that 
this requirement is indeed met, because absorption dominates over scattering by a factor of 2 to 5 at 
all densities interior to the central densities of PNSs investigated here. 

We describe neutron star matter at finite density and temperature using the Gibbs phase rules 



11]. The conditions of baryon number density and charge conservation in the mixed phase are 



n B = xnB + (1 - x) n% = X n? + (1 - x) n Q - L 



c 1 n c 1 



where % and n c are the baryon number and charge densities, respectively; H, Q, and L denote 
hadrons, quarks, and leptons. Since the dynamical time scale is much longer than the weak interaction 
time scale, beta equilibrium implies that the various chemical potentials satisfy the relations 

A* e - Hue = _ ; = biH n - qi{n e - fi Ue ) , (16) 

where hi and qi are the baryon number and electric charge of the hadron or quark of species i. When 
the neutrinos are trapped, the electron lepton number Yl = (n e + n Ue )/riB is initially fixed at a value 
~ 0.4 as suggested by collapse calculations. 

Hadronic matter is described using a field-theoretical description in which nucleons interact via the 
exchange of a, u, and p mesons. The meson-nucleon couplings and the couplings of the a self-interaction 
terms are determined by reproducing the empirical properties of nuclear matter, Eb = —16.0 MeV, 
M*/M = 0.6, K = 250 MeV, a sym = 32.5 MeV, and n = 0.16 fm~ 3 . Quark matter is described 
using the MIT Bag model, with a bag constant of B = 200 MeV/fm 3 . (Similar results are obtained 
with four-quark interactions in the Nambu-Jona-Lasinio model 0.) The mixed phase is assumed to 
be homogeneous. For more details of the calculation of the EOS, see Ref. 0. 

The cross sections per unit volume (or inverse mean free paths) for u e scattering and absorption 
are shown in Fig. ^| for the two stages of PNS evolution described before. It is important to recall 



that these curves are drawn under conditions of fixed entropy. (Constant entropy adiabats shown in 
Fig. 3 of Ref. [[J are helpful to gain insights into the behavior of the cross sections shown here.) The 
individual contributions from the different reactions in the pure nucleon and quark phases (thin lines) 
and in the mixed phase (thick lines) are marked in this figure. The vertical dashed lines show the 
central densities of 1.4M and the maximum mass configurations, respectively. Notice that quarks 
exist only in the mixed phase; the pure quark phase occurs at densities above the central densities of 
maximum mass stars in all cases shown here. 

In general, for a given density and temperature, the pure quark-phase opacity (or equivalently the 
cross section per unit volume) is less than that of hadrons due to the former's smaller matrix elements 
sampled by a relativistic phase space. In addition, for a given entropy and density, pure quark matter 
favors a lower temperature than hadronic matter ||. It is natural, therefore, that within the mixed 
phase region of hadrons and quarks, the net cross section per unit volume either flattens or decreases 
with increasing density. The reduction of opacities from that of the pure hadronic phase is enhanced 
in the v— free case, reflecting the more extreme decrease of temperature across the mixed phase region 
in that case ||. The precise density dependence of the net neutrino opacity depends upon the details 
of the mixed phase. 

Note that the total absorption cross section is larger than the scattering cross section for both 
degenerate and nondegenerate situations. As discussed above, this justifies our approximate treatment 
of scattering in the calculation of the diffusion coefficients. 

The diffusion coefficients most relevant for the PNS simulations, in matter with and without a mixed 
phase of hadrons and quarks, are D 2 and D4, and these are shown in Fig. |3|. Insofar as absorption 
dominates scattering, the behavior of these coefficients can be understood qualitatively by utilizing the 
limiting forms for neutrino absorption in the degenerate (Eq. flllD ) and nondegenerate (Eq. (|12|)) cases, 
respectively. The actual behavior is somewhat more complicated, but this assumption will suffice for 
a qualitative interpretation of Fig. |3|. In this case, the leading behaviors may be extracted to be 

D 2 oc \{fj, u /T){fj, u /T) 2 and D 4 oc X(jjl„/T = 0) , (17) 

where the mean free path A = (a/V)' 1 . In these equations, D 2 is evaluated under conditions of 
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extreme neutrino degeneracy and D4 is evaluated assuming that fi v = 0. Thus, both D2 and D4 should 
simply reflect the inverse behavior of the cross section per unit volume, which decreases with increasing 
density in the pure phases, but increases within the mixed phase region. 

Concerning the evolution of a PNS, we expect that the initial star, which is lepton rich, will not 
have an extensive mixed phase region. Only after several seconds of evolution will quark matter appear. 
In the newly-formed mixed phase region, the neutrino opacity will be substantially smaller than in the 
case in which a mixed phase region does not appear. However, due to the large v— optical depth of the 
PNS, neutrinos remain trapped, and no significant effect on emergent neutrino luminosities is expected 
at early times. As the star evolves, however, the relatively larger increase in opacity (note the increases 
in relative to D 2 ) and the growing extent of the mixed phase region eventually allows a larger flux 
of neutrinos, and thereby a more rapid evolution. 

In summary, we have calculated neutrino opacities for matter containing quarks for the tempera- 
tures, neutrino degeneracies and lepton contents relevant for PNS simulations, employing Gibbs phase 
rules to construct a mixed hadron-quark phase. We find that, in the presence of quarks, neutrinos have 
a significantly smaller opacity and hence larger diffusion coefficients than those in purely hadronic mat- 
ter at similar densities. These differences may have an observable impact on the neutrino flux from 
PNSs containing quark matter, but these differences are not expected to become apparent until the 
PNS is 10-20 s old. Simulations of PNSs with a mixed phase of hadrons and quarks are under inves- 



tigation ||16|| . The influence of heterogeneous structures [pr7| , |lq1 and superfluidity || in matter will be 
addressed in future work. 

We acknowledge research support from the U.S. Department of Energy under contracts DOE/DE- 
FG02-88ER-40388 (AWS and MP) and DOE/DE-FG02-87ER-40317 (JML). We thank Jose A. Pons 
and Sanjay Reddy for useful discussions. 
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TABLE 1: The standard model charged and neutral current vector and axial-vector couplings of 
neutrinos to leptons, baryons, and quarks; 9c is the Cabbibo angle (cos9 c = 0.973), 9w is the weak 
mixing angle (sin#yK = 0.231), and qa = 1-23 is the baryon axial-vector coupling constant. 
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FIGURE CAPTIONS 



Figure ^ u e cross sections per unit volume in pure quark matter, for T = 5 MeV. Solid lines show 



complete results from Eq. (ffl), and dashed lines indicate limiting forms from [10]. The labels "Degener- 
ate" and "Non-Degenerate" refer to neutrinos. Upper panels show results for two times nuclear matter 
density, and the lower panels show results for a neutrino energy of 50 MeV (/z is the chemical potential 
of the incoming quark). 



Figure ^ v e cross sections with various particles in matter containing a mixed phase of quarks and 
hadrons (no = 0.16 fm -3 ). The left panels show scattering cross sections for neutrinos with the indi- 
cated incoming hadrons, quarks, or leptons. Thick lines show the extent of the mixed phase region. 
The right panels show absorption cross sections on nucleons and quarks. The upper panels correspond 
to the neutrino-trapped era when s = 1 and Yl = 0.4, and the lower panels to the time following 
deleptonization when s = 2 and Y v = 0. The vertical dashed lines labelled 1*1,4 and u max indicate the 
central densities of a Mq = 1.4 M Q star and the maximum mass star {Mq = 2.22 M Q for the upper 
panels and Mq = 1.89 M for the lower panels), respectively. 

Figure |3|: Diffusion coefficients for the neutrino-trapped era (left panel) and hot deleptonized era (right 
panel). Thick lines show the extent of the mixed phase. Solid lines correspond to matter with a mixed 
phase, and dashed lines to matter containing only nucleons. The vertical dashed lines have the same 
meaning as in Fig. 2. 
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